Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

add lncRNAs annotation workflow #579

Merged
merged 12 commits into from
Mar 6, 2025
Merged

Conversation

rlibouba
Copy link
Contributor

Hello,
I would like to propose this workflow for annotating long non-coding RNAs from the GTN "Long non-coding RNAs (lncRNAs) annotation with FEELnc".

Thank you! Have a nice day!

@rlibouba rlibouba marked this pull request as draft October 22, 2024 08:32
@rlibouba rlibouba marked this pull request as ready for review November 29, 2024 10:44
Copy link
Member

@pavanvidem pavanvidem left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

It looks good but seems to follow the GTN tutorial only. The idea here is to have some generic workflows. Please make all the important parameters available to the user to set although they are not used in the tutorial.

  • make the strandness param available to set by the user
  • take both GFF3/GTF as genome annotation inputs and make the gffread step as optional if already a GTF is provided
  • Use a tiny test data

@rlibouba
Copy link
Contributor Author

Hi @pavanvidem, thank you!
I've made some changes in response to your comments.

@pavanvidem
Copy link
Member

pavanvidem commented Jan 16, 2025

thanks @rlibouba Can you please also address the following points:

  • Currently, the strandness param of the StringTie is set to "unstranded". But the input data can also be strand-specific. Please expose this parameter to the user
  • Wouldn't it possible to take both GFF3 and GTF as genome annotation inputs? StringTie can work with both types.
  • the current test data is too large (>5GB BAM file). Please subset (all the input files) to only one of the smaller scaffolds.

@rlibouba
Copy link
Contributor Author

Hello @pavanvidem,
Thanks for your feedback.
I managed to reduce the bam file a bit. I can't reduce it any more, because FEELnc doesn't have enough elements to analyze. I was able to fix the problem with the strandness parameter.

Concerning the question with StringTie, you can give it both a GFF3 and a GTF file, but this generates two outputs, one for the GFF3 and another for the GTF.

Copy link

Test Results (powered by Planemo)

Test Summary

Test State Count
Total 1
Passed 1
Error 0
Failure 0
Skipped 0
Passed Tests
  • ✅ Galaxy-Workflow-lncRNAs_annotation_workflow.ga_0

    Workflow invocation details

    • Invocation Messages

    • Steps
      • Step 1: strandedness:

        • step_state: scheduled
      • Step 2: Genome assembly:

        • step_state: scheduled
      • Step 3: Genome annotation:

        • step_state: scheduled
      • Step 4: Alignments from RNA-seq:

        • step_state: scheduled
      • Step 5: toolshed.g2.bx.psu.edu/repos/iuc/map_param_value/map_param_value/0.2.0:

        • step_state: scheduled

        • Jobs
          • Job 1:

            • Job state is ok

            Command Line:

            • cd ../; python _evaluate_expression_.py

            Exit Code:

            • 0

            Traceback:

            Job Parameters:

            • Job parameter Parameter value
              __input_ext "input"
              __workflow_invocation_uuid__ "84803768dff311efb9ab0022482354f6"
              chromInfo "/tmp/tmpdfutzrj0/galaxy-dev/tool-data/shared/ucsc/chrom/?.len"
              dbkey "?"
              input_param_type {"__current_case__": 0, "input_param": "unstranded", "mappings": [{"__index__": 0, "from": "stranded - forward", "to": "--fr"}, {"__index__": 1, "from": "stranded - reverse", "to": "--rf"}, {"__index__": 2, "from": "unstranded", "to": ""}], "type": "text"}
              output_param_type "text"
              unmapped {"__current_case__": 1, "on_unmapped": "fail"}
      • Step 6: gffread:

        • step_state: scheduled

        • Jobs
          • Job 1:

            • Job state is ok

            Command Line:

            • gffread '/tmp/tmpdfutzrj0/files/7/a/7/dataset_7a7291f4-216d-4b2a-8f72-cc0245fea1bb.dat'       -T -o output.gtf

            Exit Code:

            • 0

            Traceback:

            Job Parameters:

            • Job parameter Parameter value
              __input_ext "input"
              __workflow_invocation_uuid__ "84803768dff311efb9ab0022482354f6"
              chr_replace None
              chromInfo "/tmp/tmpdfutzrj0/galaxy-dev/tool-data/shared/ucsc/chrom/?.len"
              dbkey "?"
              decode_url false
              expose false
              filtering None
              full_gff_attribute_preservation false
              gffs {"__current_case__": 2, "gff_fmt": "gtf", "tname": ""}
              maxintron None
              merging {"__current_case__": 0, "merge_sel": "none"}
              reference_genome {"__current_case__": 0, "source": "none"}
              region {"__current_case__": 0, "region_filter": "none"}
      • Step 7: Assembly step with StringTie:

        • step_state: scheduled

        • Jobs
          • Job 1:

            • Job state is ok

            Command Line:

            • mkdir -p ./special_de_output/sample1/ &&   ln -s '/tmp/tmpdfutzrj0/files/7/a/7/dataset_7a7291f4-216d-4b2a-8f72-cc0245fea1bb.dat' guide.gff &&  stringtie '/tmp/tmpdfutzrj0/files/a/3/e/dataset_a3ee76b1-007e-4967-a7eb-91c86b0a5580.dat'  -o '/tmp/tmpdfutzrj0/job_working_directory/000/6/outputs/dataset_f91a24f4-a7a7-4948-b344-b5a82b8757f4.dat' -p "${GALAXY_SLOTS:-1}"    -G guide.gff  -b ./special_de_output/sample1/  -f '0.01' -m '200' -a '10' -j '1' -c '1' -g '50' -M '1.0'

            Exit Code:

            • 0

            Traceback:

            Job Parameters:

            • Job parameter Parameter value
              __input_ext "input"
              __workflow_invocation_uuid__ "84803768dff311efb9ab0022482354f6"
              adv {"abundance_estimation": false, "bdist": "50", "bundle_fraction": "1.0", "disable_trimming": false, "fraction": "0.01", "min_anchor_cov": "1", "min_anchor_len": "10", "min_bundle_cov": "1", "min_tlen": "200", "multi_mapping": false, "name_prefix": null, "omit_sequences": "", "point_features": null}
              chromInfo "/tmp/tmpdfutzrj0/galaxy-dev/tool-data/shared/ucsc/chrom/?.len"
              dbkey "?"
              guide {"__current_case__": 1, "coverage_file": false, "guide_source": {"__current_case__": 1, "guide_gff_select": "history", "ref_hist": {"values": [{"id": 2, "src": "hda"}]}}, "input_estimation": false, "special_outputs": {"__current_case__": 2, "special_outputs_select": "no"}, "use_guide": "yes"}
              input_options {"__current_case__": 0, "input_bam": {"values": [{"id": 3, "src": "hda"}]}, "input_mode": "short_reads"}
              rna_strandness ""
      • Step 8: lncRNA annotation step with FEELnc:

        • step_state: scheduled

        • Jobs
          • Job 1:

            • Job state is ok

            Command Line:

            • export FEELNCPATH=$(dirname $(command -v FEELnc_pipeline.sh))/../ && ln -s '/tmp/tmpdfutzrj0/files/8/1/3/dataset_813eb582-14c0-44d3-9ec5-8dc0d7b50433.dat' genome_fa && FEELnc_pipeline.sh --candidate='/tmp/tmpdfutzrj0/files/f/9/1/dataset_f91a24f4-a7a7-4948-b344-b5a82b8757f4.dat' --reference='/tmp/tmpdfutzrj0/files/1/b/5/dataset_1b554f89-8220-4974-9e44-8f574ac9d678.dat' --genome=genome_fa --outname='candidate_lncRNA' --outdir='out_feelnc'  &&  cat out_feelnc/codpot/candidate_lncRNA.codpot_RF_summary.txt out_feelnc/classifier/candidate_lncRNA.classifier.log

            Exit Code:

            • 0

            Standard Error:

            • Possible precedence issue with control flow operator at /usr/local/lib/site_perl/5.26.2/Bio/DB/IndexedBase.pm line 805.
              Filtered transcripts will be available in file: 'out_feelnc/filter/candidate_lncRNA.filter.log'
              Parsing file 'dataset_f91a24f4-a7a7-4948-b344-b5a82b8757f4.dat'...
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse input file:             [                                                                                                    ]
              Parse 
              ..
              -------                             ]
              Parse input file:             [-----------------------------------------------------------------------                             ]
              Parse input file:             [------------------------------------------------------------------------                            ]
              Parse input file:             [------------------------------------------------------------------------                            ]
              Parse input file:             [------------------------------------------------------------------------                            ]
              Parse input file:             [-------------------------------------------------------------------------                           ]
              Parse input file:             [-------------------------------------------------------------------------                           ]
              Parse input file:             [--------------------------------------------------------------------------                          ]
              Parse input file:             [--------------------------------------------------------------------------                          ]
              Parse input file:             [--------------------------------------------------------------------------                          ]
              Parse input file:             [---------------------------------------------------------------------------                         ]
              Parse input file:             [---------------------------------------------------------------------------                         ]
              Parse input file:             [----------------------------------------------------------------------------                        ]
              Parse input file:             [----------------------------------------------------------------------------                        ]
              Parse input file:             [----------------------------------------------------------------------------                        ]
              Parse input file:             [-----------------------------------------------------------------------------                       ]
              Parse input file:             [-----------------------------------------------------------------------------                       ]
              Parse input file:             [------------------------------------------------------------------------------                      ]
              Parse input file:             [------------------------------------------------------------------------------                      ]
              Parse input file:             [------------------------------------------------------------------------------                      ]
              Parse input file:             [-------------------------------------------------------------------------------                     ]
              Parse input file:             [-------------------------------------------------------------------------------                     ]
              Parse input file:             [--------------------------------------------------------------------------------                    ]
              Parse input file:             [--------------------------------------------------------------------------------                    ]
              Parse input file:             [--------------------------------------------------------------------------------                    ]
              Parse input file:             [---------------------------------------------------------------------------------                   ]
              Parse input file:             [---------------------------------------------------------------------------------                   ]
              Parse input file:             [----------------------------------------------------------------------------------                  ]
              Parse input file:             [----------------------------------------------------------------------------------                  ]
              Parse input file:             [----------------------------------------------------------------------------------                  ]
              Parse input file:             [-----------------------------------------------------------------------------------                 ]
              Parse input file:             [-----------------------------------------------------------------------------------                 ]
              Parse input file:             [------------------------------------------------------------------------------------                ]
              Parse input file:             [------------------------------------------------------------------------------------                ]
              Parse input file:             [------------------------------------------------------------------------------------                ]
              Parse input file:             [-------------------------------------------------------------------------------------               ]
              Parse input file:             [-------------------------------------------------------------------------------------               ]
              Parse input file:             [--------------------------------------------------------------------------------------              ]
              Parse input file:             [--------------------------------------------------------------------------------------              ]
              Parse input file:             [--------------------------------------------------------------------------------------              ]
              Parse input file:             [---------------------------------------------------------------------------------------             ]
              Parse input file:             [---------------------------------------------------------------------------------------             ]
              Parse input file:             [----------------------------------------------------------------------------------------            ]
              Parse input file:             [----------------------------------------------------------------------------------------            ]
              Parse input file:             [----------------------------------------------------------------------------------------            ]
              Parse input file:             [-----------------------------------------------------------------------------------------           ]
              Parse input file:             [-----------------------------------------------------------------------------------------           ]
              Parse input file:             [------------------------------------------------------------------------------------------          ]
              Parse input file:             [------------------------------------------------------------------------------------------          ]
              Parse input file:             [------------------------------------------------------------------------------------------          ]
              Parse input file:             [-------------------------------------------------------------------------------------------         ]
              Parse input file:             [-------------------------------------------------------------------------------------------         ]
              Parse input file:             [--------------------------------------------------------------------------------------------        ]
              Parse input file:             [--------------------------------------------------------------------------------------------        ]
              Parse input file:             [--------------------------------------------------------------------------------------------        ]
              Parse input file:             [---------------------------------------------------------------------------------------------       ]
              Parse input file:             [---------------------------------------------------------------------------------------------       ]
              Parse input file:             [----------------------------------------------------------------------------------------------      ]
              Parse input file:             [----------------------------------------------------------------------------------------------      ]
              Parse input file:             [----------------------------------------------------------------------------------------------      ]
              Parse input file:             [-----------------------------------------------------------------------------------------------     ]
              Parse input file:             [-----------------------------------------------------------------------------------------------     ]
              Parse input file:             [------------------------------------------------------------------------------------------------    ]
              Parse input file:             [------------------------------------------------------------------------------------------------    ]
              Parse input file:             [------------------------------------------------------------------------------------------------    ]
              Parse input file:             [-------------------------------------------------------------------------------------------------   ]
              Parse input file:             [-------------------------------------------------------------------------------------------------   ]
              Parse input file:             [--------------------------------------------------------------------------------------------------  ]
              Parse input file:             [--------------------------------------------------------------------------------------------------  ]
              Parse input file:             [--------------------------------------------------------------------------------------------------  ]
              Parse input file:             [--------------------------------------------------------------------------------------------------- ]
              Parse input file:             [--------------------------------------------------------------------------------------------------- ]
              Parse input file:             [----------------------------------------------------------------------------------------------------]
              	Your input GTF file 'candidate_lncRNA.filter.gtf' contains *113* transcripts
              	Extracting ORFs/cDNAs 1/113...
              	Extracting ORFs/cDNAs 2/113...
              	Extracting ORFs/cDNAs 3/113...
              	Extracting ORFs/cDNAs 4/113...
              	Extracting ORFs/cDNAs 5/113...
              	Extracting ORFs/cDNAs 6/113...
              	Extracting ORFs/cDNAs 7/113...
              	Extracting ORFs/cDNAs 8/113...
              	Extracting ORFs/cDNAs 9/113...
              	Extracting ORFs/cDNAs 10/113...
              	Extracting ORFs/cDNAs 11/113...
              	Extracting ORFs/cDNAs 12/113...
              	Extracting ORFs/cDNAs 13/113...
              	Extracting ORFs/cDNAs 14/113...
              	Extracting ORFs/cDNAs 15/113...
              	Extracting ORFs/cDNAs 16/113...
              	Extracting ORFs/cDNAs 17/113...
              	Extracting ORFs/cDNAs 18/113...
              	Extracting ORFs/cDNAs 19/113...
              	Extracting ORFs/cDNAs 20/113...
              	Extracting ORFs/cDNAs 21/113...
              	Extracting ORFs/cDNAs 22/113...
              	Extracting ORFs/cDNAs 23/113...
              	Extracting ORFs/cDNAs 24/113...
              	Extracting ORFs/cDNAs 25/113...
              	Extracting ORFs/cDNAs 26/113...
              	Extracting ORFs/cDNAs 27/113...
              	Extracting ORFs/cDNAs 28/113...
              	Extracting ORFs/cDNAs 29/113...
              	Extracting ORFs/cDNAs 30/113...
              	Extracting ORFs/cDNAs 31/113...
              	Extracting ORFs/cDNAs 32/113...
              	Extracting ORFs/cDNAs 33/113...
              	Extracting ORFs/cDNAs 34/113...
              	Extracting ORFs/cDNAs 35/113...
              	Extracting ORFs/cDNAs 36/113...
              	Extracting ORFs/cDNAs 37/113...
              	Extracting ORFs/cDNAs 38/113...
              	Extracting ORFs/cDNAs 39/113...
              	Extracting ORFs/cDNAs 40/113...
              	Extracting ORFs/cDNAs 41/113...
              	Extracting ORFs/cDNAs 42/113...
              	Extracting ORFs/cDNAs 43/113...
              	Extracting ORFs/cDNAs 44/113...
              	Extracting ORFs/cDNAs 45/113...
              	Extracting ORFs/cDNAs 46/113...
              	Extracting ORFs/cDNAs 47/113...
              	Extracting ORFs/cDNAs 48/113...
              	Extracting ORFs/cDNAs 49/113...
              	Extracting ORFs/cDNAs 50/113...
              	Extracting ORFs/cDNAs 51/113...
              	Extracting ORFs/cDNAs 52/113...
              	Extracting ORFs/cDNAs 53/113...
              	Extracting ORFs/cDNAs 54/113...
              	Extracting ORFs/cDNAs 55/113...
              	Extracting ORFs/cDNAs 56/113...
              	Extracting ORFs/cDNAs 57/113...
              	Extracting ORFs/cDNAs 58/113...
              	Extracting ORFs/cDNAs 59/113...
              	Extracting ORFs/cDNAs 60/113...
              	Extracting ORFs/cDNAs 61/113...
              	Extracting ORFs/cDNAs 62/113...
              	Extracting ORFs/cDNAs 63/113...
              	Extracting ORFs/cDNAs 64/113...
              	Extracting ORFs/cDNAs 65/113...
              	Extracting ORFs/cDNAs 66/113...
              	Extracting ORFs/cDNAs 67/113...
              	Extracting ORFs/cDNAs 68/113...
              	Extracting ORFs/cDNAs 69/113...
              	Extracting ORFs/cDNAs 70/113...
              	Extracting ORFs/cDNAs 71/113...
              	Extracting ORFs/cDNAs 72/113...
              	Extracting ORFs/cDNAs 73/113...
              	Extracting ORFs/cDNAs 74/113...
              	Extracting ORFs/cDNAs 75/113...
              	Extracting ORFs/cDNAs 76/113...
              	Extracting ORFs/cDNAs 77/113...
              	Extracting ORFs/cDNAs 78/113...
              	Extracting ORFs/cDNAs 79/113...
              	Extracting ORFs/cDNAs 80/113...
              	Extracting ORFs/cDNAs 81/113...
              	Extracting ORFs/cDNAs 82/113...
              	Extracting ORFs/cDNAs 83/113...
              	Extracting ORFs/cDNAs 84/113...
              	Extracting ORFs/cDNAs 85/113...
              	Extracting ORFs/cDNAs 86/113...
              	Extracting ORFs/cDNAs 87/113...
              	Extracting ORFs/cDNAs 88/113...
              	Extracting ORFs/cDNAs 89/113...
              	Extracting ORFs/cDNAs 90/113...
              	Extracting ORFs/cDNAs 91/113...
              	Extracting ORFs/cDNAs 92/113...
              	Extracting ORFs/cDNAs 93/113...
              	Extracting ORFs/cDNAs 94/113...
              	Extracting ORFs/cDNAs 95/113...
              	Extracting ORFs/cDNAs 96/113...
              	Extracting ORFs/cDNAs 97/113...
              	Extracting ORFs/cDNAs 98/113...
              	Extracting ORFs/cDNAs 99/113...
              	Extracting ORFs/cDNAs 100/113...
              	Extracting ORFs/cDNAs 101/113...
              	Extracting ORFs/cDNAs 102/113...
              	Extracting ORFs/cDNAs 103/113...
              	Extracting ORFs/cDNAs 104/113...
              	Extracting ORFs/cDNAs 105/113...
              	Extracting ORFs/cDNAs 106/113...
              	Extracted '106' ORF/cDNAs sequences on '113'.
              > Run random Forest on '/tmp//413_candidate_lncRNA.codpot.test_rna.fa'
              	1. Compute the size of each sequence and ORF
              	2. Compute the kmer ratio for each kmer and put the output file name in a list
              	3. Compute the kmer score for each kmer size on learning and test ORF
              	4. Merge the score and size files into one file for each type
              	5. Make the model on learning sequences and apply it on test sequences
              		Running 10-fold cross-validation on learning
              
                |                                                                            
                |                                                                      |   0%
                |                                                                            
                |=======                                                               |  10%
                |                                                                            
                |==============                                                        |  20%
                |                                                                            
                |=====================                                                 |  30%
                |                                                                            
                |============================                                          |  40%
                |                                                                            
                |===================================                                   |  50%
                |                                                                            
                |==========================================                            |  60%
                |                                                                            
                |=================================================                     |  70%
                |                                                                            
                |========================================================              |  80%
                |                                                                            
                |===============================================================       |  90%
                |                                                                            
                |======================================================================| 100%
              		10-fold cross-validation step is finish. Best threshold found: '0.6356'.
              		Write the coding label for '413_candidate_lncRNA.codpot.modelTest.out' in 'out_feelnc/codpot//candidate_lncRNA.codpot_RF.txt'.
              > Writing the GTF output files
              
              Transcripts without ORF found: check file 'out_feelnc/codpot//candidate_lncRNA.codpot.noORF.gtf'.
              	These transcripts correspond probably to lncRNAs since no ORF was found. You might also want to change the --testorftype option.
              window : 10000 - max window : 100000 - lncrna : out_feelnc/codpot//candidate_lncRNA.codpot.lncRNA.gtf - mrna : /tmp/tmpdfutzrj0/files/1/b/5/dataset_1b554f89-8220-4974-9e44-8f574ac9d678.dat - biotype: 0
              Possible precedence issue with control flow operator at /usr/local/lib/site_perl/5.26.2/Bio/DB/IndexedBase.pm line 805.
               Importation  merge GTF file : 106 lncRNA done 
               Importation  merge GTF file : 13795 mRNA done 
               ... Looking for interactions, currently eq to 300 
              

            Standard Output:

            • ####################
              Create output directory out_feelnc
              
              ####################
              Step 1/3: Running the FEELnc filter module:
              
              commande line:
              	FEELnc_filter.pl --infile=/tmp/tmpdfutzrj0/files/f/9/1/dataset_f91a24f4-a7a7-4948-b344-b5a82b8757f4.dat --mRNAfile=/tmp/tmpdfutzrj0/files/1/b/5/dataset_1b554f89-8220-4974-9e44-8f574ac9d678.dat --biotype=transcript_biotype=protein_coding --monoex=-1 --outlog=out_feelnc/filter/candidate_lncRNA.filter.log > out_feelnc/filter/candidate_lncRNA.filter.gtf
              
              
              ####################
              Step 2/3: Running the FEELnc codpot module:
              
              commande line:
              	FEELnc_codpot.pl --infile=out_feelnc/filter/candidate_lncRNA.filter.gtf --mRNAfile=/tmp/tmpdfutzrj0/files/1/b/5/dataset_1b554f89-8220-4974-9e44-8f574ac9d678.dat --genome=genome_fa -b transcript_biotype=protein_coding --mode=shuffle --outname=candidate_lncRNA.codpot --outdir=out_feelnc/codpot/
              
              
              ####################
              Step 3/3: Running the FEELnc classifier module:
              
              commande line:
              	FEELnc_classifier.pl --lncrna=out_feelnc/codpot//candidate_lncRNA.codpot.lncRNA.gtf --mrna=/tmp/tmpdfutzrj0/files/1/b/5/dataset_1b554f89-8220-4974-9e44-8f574ac9d678.dat --log=out_feelnc/classifier/candidate_lncRNA.classifier.log > out_feelnc/classifier/candidate_lncRNA.classifier.txt
              
              
              
              ####################
              Done: results in out_feelnc/{filter;codpot;classifier}
              
              # Summary file:
              -With_cutoff:	0.6356
              -Nb_lncRNAs:	106
              -Nb_mRNAs:	NA
              #FEELnc Classification
              #lncRNA file :  lncrna : out_feelnc/codpot//candidate_lncRNA.codpot.lncRNA.gtf 
              #mRNA file : /tmp/tmpdfutzrj0/files/1/b/5/dataset_1b554f89-8220-4974-9e44-8f574ac9d678.dat
              #Minimal window size : 10000
              #Maximal window size : 100000
              #Number of lncRNA : 106 
              #Number of mRNA : 13795
              #Number of interaction : 695 
              #Number of lncRNA without interaction : 0
              #List of lncRNA without interaction : 
              

            Traceback:

            Job Parameters:

            • Job parameter Parameter value
              __input_ext "input"
              __workflow_invocation_uuid__ "84803768dff311efb9ab0022482354f6"
              chromInfo "/tmp/tmpdfutzrj0/galaxy-dev/tool-data/shared/ucsc/chrom/?.len"
              dbkey "?"
      • Step 9: Concatenation step:

        • step_state: scheduled

        • Jobs
          • Job 1:

            • Job state is ok

            Command Line:

            • python /tmp/tmpdfutzrj0/galaxy-dev/tools/filters/catWrapper.py '/tmp/tmpdfutzrj0/job_working_directory/000/8/outputs/dataset_263acba7-c500-46ef-82e6-ddd2e8c74464.dat' '/tmp/tmpdfutzrj0/files/e/d/a/dataset_eda56677-43d0-43a3-bcbd-51e9507a3f2d.dat' '/tmp/tmpdfutzrj0/files/1/b/5/dataset_1b554f89-8220-4974-9e44-8f574ac9d678.dat'

            Exit Code:

            • 0

            Traceback:

            Job Parameters:

            • Job parameter Parameter value
              __input_ext "gtf"
              __workflow_invocation_uuid__ "84803768dff311efb9ab0022482354f6"
              chromInfo "/tmp/tmpdfutzrj0/galaxy-dev/tool-data/shared/ucsc/chrom/?.len"
              dbkey "?"
              queries [{"__index__": 0, "input2": {"values": [{"id": 5, "src": "hda"}]}}]
    • Other invocation details
      • history_id

        • 5982786cc98d0e7f
      • history_state

        • ok
      • invocation_id

        • 5982786cc98d0e7f
      • invocation_state

        • scheduled
      • workflow_id

        • 5982786cc98d0e7f

@pavanvidem
Copy link
Member

hi @rlibouba thanks for adapting the workflow.

  • please add a release
  • make sure that all the inputs have the annotation
  • please check that no underscrores are allowed in the inputs and output names

Copy link
Member

@mvdbeek mvdbeek left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thanks so much @rlibouba!

@mvdbeek mvdbeek merged commit 3a1a9d0 into galaxyproject:main Mar 6, 2025
7 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants